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To  understand  and  control  the  basic  energy  transport 
phenomena  in  the  Vapor  Shield  at  different  material 
surfaces  under  high  heat  flux. 

To  explore  the  effect  of  a  strong  magnetic  field  in 
decreasing  surface  erosion  and  drag  coefficients  under 
high  heat  flux  under  EM  gun  conditions. 

To  model  the  Magnetic  Vapor  Shield  mechanism  and  to 
compare  with  experimental  results. 

To  understand,  characterize,  and  control  the  plasma- 
propellant  energy  transport  process  by  measuring 
propellant  bum  rates. 

To  explore  and  model  the  plasma-propellant  interaction 
process  and  the  flame  vapor  shield  mechanism. 

To  systematically  diagnose  and  explore  the  reaction  of 
different  propellants  under  different  injected  plasma 
conditions. 


Majoir  Issii©s  ft®  Imv©sfti§aft© 

□  To  systematically  explore  the  resistance 
of  various  materials  (insulators, 
composites,  metals  and  alloys,  non- 
metallic  conductors,  metallic  and 
diamond  coatings,  etc.)  to  high  heat 
flux  erosion,  and  plasma  reactions  with 
solid  and  liquid  propellants. 

□  To  determine  the  nature  of  drag  forces 
and  turbulent  convective  flow  in 
contributing  to  surface  erosion. 

□  To  fully  diagnose  the  vapor  shield 
plasma,  and  the  plasma  and  heat  flux 
distribution  onto  propellant  samples. 

□  To  model  and  predict  the  vapor  shield 
effectiveness  both  with  and  without  a 
magnetic  field,  and  model  and  predict 
the  flame  vapor  shield  effect  and  burn 
rate  of  propellants. 


□  Predicted  vapor  shield  transmission  factors  that  were 
verified  by  experiment.  Energy  transmission  factor 
through  the  vapor  shield  is  on  the  order  of  10%. 
Radiation  transport  is  the  most  important  energy  transfer 
mechanism.  Demonstrated  that  some  materials  have 
better  vapor  shielding  properties  than  others. 

□  Erosion  of  different  materials  was  explored,  and  the 
axial  erosion  dependence  was  measured  for  different 
materials.  SEM,  EDXA  and  Auger  analyses 
demonstrate  the  redeposition  of  electrode  materials  at 
different  locations.  Conditioning  effects  of  single 
material  samples  exposed  to  multiple  shots  were  noted, 
which  will  be  important  in  operation  of  railguns. 

□  The  spectra  of  observed  plasma  light  emission  indicated 
the  presence  of  electrode  and  sample  materials. 

□  Preliminary  results  on  the  effect  of  the  magnetic  field. 
Magnetic  fields  decreased  surface  erosion  by  20%  for 
fields  of  5  Tesla. 

□  Developed  series  of  1-D,  time-dependent,  MHD  codes  to 
model  vapor  shield  phenomena,  ZEUS,  a  global,  time 
dependent  code  for  the  source  plasma,  and  MAGFIRE,  a 
1-D,  time  dependent  MHD  code  with  coupled  radiation 
transport  for  ablative  plasma. 
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□  Generated  materials  erosion  data  base  for  metals,  alloys, 
insulators,  non-metallic  conductors  and  refractory 
materials.  Demonstrated  that  metallic  surfaces  have 
mixed  melting/erosion  with  strong  axial  erosion 
dependence,  while  insulators  and  graphites  ablate 
homogeneously  along  the  axis. 

□  Modified  the  chamber  geometry  (10  inch  cross  connected 
to  the  6  inch  cross). 

□  Modified  diagnostics  arrangements  and  developed  in¬ 
bore  heat  flux  (direct  and  indirect  measurement  of  heat 
flux  to  the  wall)  with  1-D  and  2-D,  time  dependent  heat 
flux  codes  “SURFHEAT  and  2DSURFHEAT’’ ,  and 
plasma  conductivity  diagnostic  technique. 

□  Obtained  empirical  general  scaling  law  of  the  energy 

transmission  factor  /  through  the  vapor  shield  layer  (for 
different  materials)  for  ideal  and  non-ideal  plasma 
models,  /  =  Ki  +  K2  where  Ki  and  K2  are 

constants  depending  on  material  properties,  S  is  the 
source  fluence. 

□  Obtained  results  on  the  effect  of  the  magnetic  field 
(Magnetic  Vapor  Shield)  with  fields  up  to  8.75T. 
Magnetic  fields  decreased  surface  erosion  by  up  to  35% 
for  fields  of  6.25  Tesla.  An  empirical  general  scaling 
law  of  the  energy  transmission  factor  for  Lexan  with  the 
effect  of  the  magnetic  field,  /(s,B)  has  been  obtained, 

/(S,B)  =  ^l(B)  S-1/4+  <E>2(B),  where  ^l(B),  <h2(B)  are 
quadratic  functions  of  the  magnetic  field. 


□  Identified  species  that  appear  in  the  majority  of  the 
discharges.  C2  (Swan  bands),  CH4  (absorption  line), 
neutral  copper  (from  the  electrodes),  and  neutral  carbon. 

□  Measured  the  source  average  plasma  temperature  using 
optical  emission  spectroscopy,  correlated  to  conductivity 
probes  and  heat  flux  measurements,  and  calculations 
based  on  the  source  ablation  and  ZEUS -code  predictions. 

□  Continuous  materials  erosion  studies  on  coated  surfaces 
with  different  coating  techniques,  tungsten,  tantalum  and 
chromium  alloys,  ceramic  insulators  and  diamond 
coatings. 

□  Modified  SURFHEAT  code  to  include  melting.  Plasma 
temperature  deduced  from  heat  flux  measurements 
correlates  well  with  values  obtained  from  optical 
emission  spectroscopy.  The  updated  version  of 
SURFHEAT  yields  more  accurate  results  in  heat  flux 
diagnostics. 

□  Detailed  analysis  of  the  magnetic  vapor  shield,  including 
convection  heat  flux,  velocity  slowdown,  and  photon 
attenuation  within  the  boundary  layer. 

□  Accelerating  projectiles,  0.33  g  Lexan  projectiles,  and 
velocity  measurements  using  break  wires,  photo¬ 
interrupters,  laser  cut-off,  and  absolute  pressure 
transducers  “fuse-atmospheric  and  fuseless-vacuum 
operation”. 

□  Drag  forces  measurements,  using  absolute  pressure 
transducers  inside  a  special  diagnostics  barrel. 
Measured  0.8  kbar  pressure  inside  the  barrel  at  5  kJ 
input  energy  “fuse  and  fuseless  operation”. 


DISCHARGE  VOLTAGE  (10  kV  max.)  1  -  8  kV 


PEAK  CURRENT  (100  kA  max.) 

NET  INPUT  ENERGY  (15  kJ  max.) 

DISCHARGE  PERIOD 

RADIATED  POWER 

PEAK  PRESSURE 

PLASMA  DENSITY 

PEAK  PLASMA  TEMPERATURE 

AVERAGE  PLASMA  TEMPERATURE 


20-100kA 
1  -  8  kJ 
100  4sec 
2-70  GW/m2 
>  1  kbar 
lO^^-lO^^m"^ 
4  -  6  eV 
1  -  3  eV 


AVERAGE  PLASMA  VELOCITY 


«12  km/sec 


Diagnostic  Ports 


Inches 


2400^F 

I  II  III 


TOP;  Schematic  of  SIRENS  showing  the  source  section 

(capillary),  barrel  section  with  magnet  assembly 
(magnetic  vapor  shield  studies). 

BOTTOM:  Diagnostics  arrangement  showing  standard  diagnostics 
(currents,  voltage,  B-dots),  fiber  optics  for  optical 
emission  spectroscopy,  and  a  thermocouple  for  heat  flux 
measurements. 


300|iF 

SOURCE  SECTION  BARREL  SECTION  TARGET 


Optomteirapters 


SOURCE  SECTION  BARREL  SECTION  TARGET 

ASSEMBLY 

BARREL  DIAGNOSTICS 


TOP:  conductivity  probes  (resistivity),  thermocouples  (heat 

flux),  and  pressure. 

BOTTOM:  Break  wires  (velocity),  pressure  (drags),  and 
optointerrupters  (velocity). 


Photographs  of  SIRENS  showing  the  source  input 
section  and  magnet  current  feedthrough  (top  photo),  and 
target  and  expansion  chamber  with  diagnostics 
feedthrough  (bottom  photo). 
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Discharge  Current  (kA) 


Shot  #  362 

Net  Energy  Input:  5.008  kJ 
Material:  Tantalum  on  Copper  Substrate 
(Ta  sputtering  on  Cu  substrate,  25  p  thick) 
“RUN  No.  125-4" 


0  10  20  30  40  50  60  70  80  90  100 


Discharge  time  (jisec) 


Oscillograms  of  discharge  current  and  voltage  (top),  and  calculated  power  and 

resistance  (bottom). 


Plasma  Resistance  (Q) 


Temperature  (Kelvin) 


INCIDENT 

PLASMA 


SCHEMATIC  OF  HEAT  FLUX  MEASUREMENTS 
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MEASURED  AND  CALCULATED 
TEMPERATURE  HISTORIES 


2  T*  T* 

_ i2 


At 


withRjp  = 


kA 


jJ 


=  Thennal  resistance  between  nodes  j  and  p 


]p 


Nodal  equations: 


Front  surface :  t‘"‘  =  t,‘ + -  ,{k(T‘  -  t,‘J  +  ax  q'{ 

pC„  Ax 


Interior  region:  C'=  Tl+~p^lk{T:.,  -  +  -T;jf 


Back  surface:  t;:;»  n.,+j|^lk(T;., -i^J+Axq' 


Time  (sec) 


COMPARISON  BETWEEN  TEMPERATURE 
HISTORIES  OBTAINED  BY  SLOW  AND  FAST 
RESPONSE  THERMOCOUPLES 
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Input  Energy  (kJ) 


^  Ablation  Algorithm 
□  Tungsten  Target 
O  Pyrolytic  Graphite  Target 
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HEAT  FLUX  RESULTS  FOR  A  TUNGSTEN  ALLOY 
TARGET  (MELT  LAYER  REMOVED)  COMPARED 
TO  PYROLYTIC  GRAPHITE  TARGET  (ABLATION 
OF  TARGET  SURFACE),  AND  SOURCE  ABLATION 
ALGORITHM. 


t  (^sec) 

Conductivity  probes  are  connected  to  a  modular  potential  divider,  where  the 
potential  between  each  two  probes  is  recorded.  Probe  potentials,  discharge 
current  and  voltage  data  are  recorded  with  a  LeCroy  6810  waveform  digitizer 
interfaced  to  a  PC  via  GPIB/PC-2.  Through  data  analysis  and  manipulation, 
plasma  resistivity  (and  consequently  plasma  temperature)  can  be  extracted. 


PLASMA  RESISTANCE  OBTAINED  FROM  CONDUCTIVITY  PROBES 
ALONG  THE  BARREL  AXIAL  DIRECTION 


The  resistivity  model  for  weakly  nonideal  plasmas  is 
evaluated  from: 

ln(l  +  1.4A^)l/2  ^ 

ei  eiO  ln(A)  -  r 

where.  In  (  1 +  1.4  replaces  the  Coulomb 

logarithm  ln(A)  to  account  for  the  nonideal  effects,  and 
Am  =  (A  Ts*)  where  rs*(=  ts/Xd)  is  the  normalized 
non-Debye  screening  radius;  ^eiO  is  the  Spitzer 
resistivity  model,  and  €1  represents  the  term  for  the 
localization  of  electrons  diue  to  incomplete  screening: 

a  =  ( 1  +  0.8  —  )^ 

where  w  is  the  energy  dispersion  due  to  electric 
micro  fields. 


Plasma  Resistivii? 


Ideal  and  nonideal  plasma  resistivities  for  a  plasma 
density  of  and  an  average  charge  state  =1.0 


-LLl  J._l  1-4  4  t  l 


Plasma  temperature  as  a  function  of  the  discharge  time,  calculated  from 
conductivity  probes  (CPl-2  and  CPl-3)  and  discharge  current,  for  1.17  kJ 
input  energy. 


Plasma  Tempeiatnie  (eV) 
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2  3  4 

Input  Energy  (kJ) 

Temperatures  of  the  source  plasma  deduced  from  ablation 
measurements  and  ZEUS  code.  Also  shown  are  temperatures  as 
calculated  by  Optica]  Emission!  Spectroscopy 


Input  Energy  (kJ) 


Average  plasma  temperatures  in  the  barrel  deduced  from  heat  flux 
and  conductivity  measurements.  Also  shown  temperatures  as 
determined  by  Optical  Emission  Spectroscopy. 
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Comparison  between  measured  erosion  for  pure  metals,  alloys,  coated 
surfaces,  insulators  and  graphite,  at  an  incident  heat  flux  of  33  GW/m2  over 


Erosion 


I LC  Chromium  i 
1  Electroplating 

LC  Chromium! 
Sputtering 

HC  Chromium 
Electroplating 

HC  Chromium 
Sputtering 

Comparison  of  L.C.  and  H.C.  Chromium  coatings  on 
copper,  for  tested  samples  prepared  via  electroplating 
and  sputter  deposition. 
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Molybdenum/Copper 


Tantalum/Copper 


■  Calculated  (no  heat  sink) 

■  Calculated  (with  heat  sink) 

■  Measured 


j 


Calculated  (with  and  without  heat  sink)  and  measured 
erosion  depths  of  molybdenum  and  tantalum  on  copper 
substrates  exposed  to  33  GW/m2  heat  flux.  Copper 
substrates  have  6.35  mm  diameter  and  3.175  mm 
thickness.  Molybdenum  and  tantalum  coatings 
(prepared  by  sputtering)  are,  approximately,  of  30  pm 
thickness.  Calculated  values  are  based  on  full  energy 
deposition  without  vapor  or  melt  shield. 


■  Calculated  (No  HS,  No  VS) 
S  Measured 

■  Calculated  with  HS 


Calculated  and  measured  erosion  of  pure  tungsten,  tungsten  coating  on  molybdenum  and  copper, 
and  molybdenum  and  tantalum  coating  on  copper.  CVD  =  Chemical  Vapor  Deposition, 
LPPS  =  Low  Pressure  Plasma  Spray,  MSD  =  Molten  Salt  Deposition.  HS  and  VS  stand  for  Heat 
Sink  and  Vapor  Shield,  respectively. 


heat  fluence  (SIRENS  experiment). 


Oat»s/Otx>t. 


100 


VAPOR  SHIELD  EFFECT 

Comparison  between  SIRENS  and  Efrimov 
results  on  graphite 


Energy  Transmission  factor 
through  the  vapor  shield 


HF  GW/m2 


Energy  transmission  factor  through  the 
boundary  layer  “vapor  shield”  for  three 
insulators  of  interest  for  EM,  ET  and  ETC 
launchers 


Individual  and  combined  effects  of  plasma  velocity 
slowdown,  attenuated  photon  intensity  and  suppression 
of  turbulence  on  surface  ablation  as  functions  of  the 
externally  applied  parallel  magnetic  field,  at  an  input 
energy  of  3  kJ.  The  combined  effect  on  ablation  is  in 
good  agreement  with  the  experimentally  obtained 
normalized  ablation. 
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Input  Energy  (kJ) 

Peak  pressures  for  fuseless  (vacuum) 
and  fuse  operation  (atmospheric) 
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•  P2  max 
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Wire  Circuit  Voltage  (mV) 


SHOT  #  494 
INPUT  ENERGY  1.2  KJ 
FUSE  OPERATION  “ATMOSPHERIC” 
BREAK-WIRES  TRACES 


Free  Expansion 
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Fuse  Operation  (Atmospheric) 

Projectile:  330  mg  Lexan,  6  mm  diameter 


Shot  #  0414,  Sample  #  W220 
Incident  Heat  Flux  33  GW/m2 


Unexposed  sample,  magnification  2000x 


Exposed  sample,  magnification  3000x 


SEM  micrographs  of  unexposed  (top)  and  exposed 
(bottom)  diamond  coating  on  a  single  crystal  silicon 
wafer  (10-13  ^m  diamond  coating  thickness  on  25.4  mm 
diameter  x  0.5  mm  thick  single  crystal  silicon  wafer, 
coating  processing  parameters:  2%  CH4,  55  Torr,  300 
scans,  26  hours). 


Shot  #  0415,  Sample  #  W26 
Incident  Heat  Flux  33  GW/m2 


Exposed  sample,  magnification  2000x 


SEM  micrographs  of  unexposed  (top)  and  exposed 
(bottom)  diamond  coating  on  a  single  crystal  silicon 
wafer  (10-13  ^im  diamond  coating  thickness  on  25.4  mm 
diameter  x  0.5  mm  thick  single  crystal  silicon  wafer, 
coating  processing  parameters:  2%  CH4,  55  Torr,  300 
scans,  26  hours). 


Shot  #  0464,  Sample  #  PYl 
Incident  Heat  Flux  4  GW/m^ 


Uaexposed  sample,  magnification  500bx 


Exposed  sample,  center  of  impact,  magnification  5000x 


SEM  micrographs  of  unexposed  (top)  and  exposed 
(bottom)  diamond  coating  on  a  single  crystal  silicon 
wafer  (15  pm  diamond  coating  thickness  on  25.4  mm 
diameter  x  3.3  mm  thick  single  crystal  silicon  wafer, 
coating  processing  parameters:  2%  CH4,  55  Torr,  300 
scans,  32  hours). 


Shot  #  0470,  Sample  #  W-llN 
Incident  Heat  Flux  3  GW/m2 


Unexposed  sample,  magnification  3000x 
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Exposed  sample,  magnification  3000x 


SEM  micrographs  of  unexposed  (top)  and  exposed 
(bottom)  diamond  coating  on  a  single  crystal  silicon 
wafer  (20  ^m  diamond  coating  thickness  on  25.4  mm 
diameter  x  3.3  mm  thick  single  crystal  silicon  wafer, 
coating  processing  parameters:  2%  CH4,  55  Torr,  300 
scans,  36  hours). 


Shot  #  0471,  Sample  #  W-12N 
Incident  Heat  Flux  1  GW/m2 


Unexposed  sample,  magnification  3000x 


Exposed  sample,  edge  of  impact,  magnification  3000x 


SEM  micrographs  of  unexposed  (top)  and  exposed 
(bottom)  diamond  coating  on  a  single  crystal  silicon 
wafer  (20  pm  diamond  coating  thickness  on  25.4  mm 
diameter  x  3.3  mm  thick  single  crystal  silicon  wafer, 
coating  processing  parameters:  2%  CH4,  55  Torr,  300 
scans,  36  hours). 


^  Continuous  materials  erosion  studies,  potentially  on  full 
section  of  a  railgun  barrel. 


Acceleration  studies  using  several  projectile  geometry 

and  material,  on  fuseless  “vacuum”  and  fuse 
“atmospheric”  operation. 


Drag  forces  studies  with  and  without  a  projectile. 


Development  of  new  diagnostics  for  boundary  layer  and 
armature  physics  studies. 


®  Studies  on  plasma-propellant  interactions  “solid  and 
liquid  propellants”  and  energy  transport  process 
“propellant  burn  rates”,  and  the  flame  vapor  shield 
mechanism. 

^  Modeling  and  prediction  of  the  flame  vapor  shield  effect 
and  burn  rate  of  propellants. 


^  Detailed  measurements  and  analysis  of  the  magnetic 

vapor  shield  effect  at  higher  values  of  magnetic  fields, 
and  higher  heat  fluxes. 

^  Ablation  physics  at  surface  “boundary  layer”  (1-D), 
includes  heat  transfer. 


^  Turbulent  plasma  boundary  layer  analysis  (2-D),  to 
explain  magnetic  field  effects. 


Hatezial 
Test  Stand 


ISOMETRIC  AND  CUTAWAY  VIEWS  OF  THE  MODIFIED 
COMBUSTION  CHAMBER,  SHOWING  DIAGNOSTICS  ACCESS, 
PROPELLANT  AND  SAMPLE  TEST  STANDS  (PIVOTING),  AN 
ABSOLUTE  PRESSURE  TRANSDUCER  FOR  CHAMBER 
PRESSURE,  AND  THE  LOADING  AND  SERVICES  PORT  WITH  A 
PRESSURE  RELIEF  VALVE. 


Ycrtical 

IniectioiL 

Source 


0  12  3  4 

(Incites) 


CONCEPTUAL  DESIGN  OF  THE  MODIFIED  COMBUSTION 
CHAMBER.  THE  CHAMBER  IS  CUBICAL  “6  INCH,  6-WAY 
STAINLESS  STEEL  CUBE”  AND  INCLUDES  TWO  3-D 
POSITIONING  PIVOT  TEST  STANDS  FOR  THE  PROPELLANT 
AND  THE  MATERIAL  SAMPLES. 


Glidcop 
rail 


